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Abstract 

We argue that as the atomic number of the target nucleus A — > oo, the mul- 
tiplicity of leading particles in hadron- nucleus collisions tends to a finite limit. 
The limiting multiplicities for various particle production are computed for 
both proton and pion projectiles. Signatures at finite A are discussed. Data 
from 100 GeV/c central hadron- nucleus collisions are analyzed and found to 
be in qualitative agreement with this picture. 
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The production of leading hadrons in high-energy hadron-nucleus scattering has been 
studied for many years (for a recent review, see |I|]). The inclusive data indicate a strong 
suppression of leading particles in central hadron-nucleus (h-A) collisions. Even for heavy 
nuclei, the main contribution to the inclusive spectra comes from the interactions near the 
nuclear edge and not the center. More recent attempts to enhance the contribution of 
collisions at small impact parameters by selecting events with multiple nucleon emission 
confirm that the leading hadron multiplicity rapidly decreases with increasing A. 

There are few theoretical results for the leading particle multiplicity from central h-A 
collisions. The earliest phenomenological results were based on the applications of various 
models of soft strong interactions, notably Gribov reggeon calculus. There are some clear 
problems with applying this approach to the projectile fragmentation region (which we will 
discuss below) and in fact very few explicit treatments of the leading particle multiplicity 
problem have been made. An alternative approach is to apply perturbative QCD (pQCD), 
in which parton correlators are the elementary building blocks, see e.g. 0. There have 
been theoretical arguments |Q, and more recently experimental evidence from nuclear size 
dependence tests , which show that hadronization of high-energy partons occurs outside 
the nucleus. 

The new observation of this paper is that for large nuclei, and a sufficiently high energy 
projectile, the origin of leading particles may have a very simple explanation based on the 
underlying parton dynamics and may be calculable through the pQCD ingredients of the 
strong interaction dynamics. 

Before addressing our proposed pQCD mechanism, let us first review the expectations 
of the Gribov-Glauber theory of hadron-nucleus interactions 0. This theory is known to 
describe quantitatively the total and elastic cross sections of high-energy hadron-nucleus 
scattering, see [0] and references therein. Corresponding reggeon calculus diagrams include 
multipomeron exchanges with a hadron, which are expressed through the vertices where 
the hadron couples to n pomerons. These vertices involve complicated hadron intermedi- 
ate states and hence are not known independently. The use of the Abramovskii, Gribov, 
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Kancheli (AGK) cutting rules allows one to describe also the production of hadrons in 
the central rapidity range, for a review see e.g. |^. However, the production of particles in 
the projectile fragmentation region requires additional information because the AGK cutting 
rules are not applicable for this kinematic region |^. A generic Gribov-Glauber approach 
does not have predictive power here. One has to know the structure of the vertex (so called 
Mandelstam cross) for coupling of the hadron with n pomerons, including information on 
how the energy of the projectile is split between several pomerons. One example of such 
a model is an eikonal type model [l^ where it is assumed that the energy is split equally 
between n interactions. Another is the additive quark model where the eikonal is applied to 
the interactions of constituent quarks [jlT|. One common feature of these models is that in 
the limit of A — ^ oo the leading particle multiplicity tends to zero. 

At the same time according to QCD-parton model concepts, leading particles originate 
from the emerging fast partons of the collision debris. Here a leading or fast parton is defined 
as one carrying a large fraction of the projectile's longitudinal momentum. In particular 
there is a wide rapidity separation between a fast parton and the sea partons. A reasonable 
criterion is that x > 0.2, where x is its momentum fraction. In addition, the absolute 
momentum of the fast parton should be sufficiently large so that it cannot easily interact 
softly with the medium. According to the space-time picture of strong interactions ||12| a 



parton fluctuates into other hadronic states in a time governed by the uncertainty principle. 
Thus a parton of three-momentum p fluctuates into a state of mass m in a time, 

(1) 

The state of lowest mass gives the characteristic time required for the point-like parton to 
become spatially extent. Conservatively taking m ~ rup (versus m^r) one finds that the 
condition, 

t=^>2RA (2) 

is satisfied at A=200 for p > 15 GeV/c, where Ra is the radius of a nucleus of atomic 
number A. This requires projectiles with energy E > 75 GeV. 



In the ideal parton model limit, such fast partons would interact rarely with the sur- 
rounding nuclear medium. In this ideal description, the interaction of the hadron projectile 
with the nuclear target would be primarily through the wee partons in the former. As such, 
a fast parton or a coherent configuration of fast partons would filter through essentially 
unaltered. This would imply that the leading particle spectrum for a given hadronic pro- 
jectile on any nuclear target would be the same Such universality in the spectra is 
qualitatively inconsistent with experiment 1^. 

QCD also predicts that fast partons will lose a negligible fraction of their longitudinal 
momentum if pinc — > oo and A = const. This expectation for fast partons to retain their 
longitudinal momentum can be tested from the nuclear size dependence of the Drell-Yan 
cross section. Recall that Drell-Yan production depends kinematically only on the longi- 
tudinal momentum fractions of the two impinging partons. Thus if there is no substantial 
attenuation of longitudinal momentum for a fast parton, the Drell-Yan cross section will 



scale linearly with A. This is what experiment finds ||14|| . 

In addition, QCD predicts [|T^ that fast partons undergo quasi-elastic rescatterings which 
give them a transverse momentum < > that increases linearly with the distance traveled 
in the nuclear medium. Up to geometrical corrections the dependence goes approximately 
as A^^^. This QCD prediction has been confirmed by the Drell-Yan experiment ||14|. The 



Pt - broadening in the dimuon spectrum for hadron projectile hp and target nucleus A, 
^{PT)hpA, is given by the transverse momentum of the initiating partons in the projectile. 
The effect of coherence between collisions of the fast parton with the nuclear media does 
not change the hnear dependence on distance for the average {pT)hpA, as discussed recently 
in |6|. 



In pQCD transverse broadening and energy loss are related as [|T^[] 

AE ^ ^A(p|)L/2, (3) 



^This relation depends weakly on the approximations used in ||T^, see discussion in [17|. 
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where L is the path in the nuclear matter. Transverse momentum smearing imphes that 
a type of hmiting curve is expected. This hmiting behavior arises because an increase in 
the relative transverse momentum between two leading partons decreases the probability 
for two such partons to coalesce. In particular one expects a high coalescence probability 



for relative transverse momenta = y'^ipr) that are of typical hadronic scale, e.g. pt^ < 
^typical ^ y^2(0.3)^ GcV/c ~ 0.420 GeV/c. Increasing the nuclear volume, thus increasing 
PTr,, decreases the likelihood of coalescence at least as oc ^ for large pt, see equation (8) 
below. So in the limit of infinite nuclear volume, A oo, the coalescence probability goes to 
zero, and the leading partons should fragment independently and so produce independent 
jets. At the same time the fractional energy loss tends to zero for a fixed large value of 
A. Thus the ^-distribution of these jets will not depend on A, although the transverse 
momentum will increase with A. Obviously one can chose the double limit of E oc A'', 
A ^ oo so that p'^ > cA"^ while both condition (Q) and the inequality ^ — > are satisfied. 
This corresponds to 

k>n + 1/3. (4) 

In this limit, when the leading partons (mostly quarks in a nucleon projectile or q and q 
in a pion) fragment independently, it is possible to calculate the leading parton production 
cross section integrated over the transverse momentum pt- Let zda^'^^ {z) / dz denote the 
one-particle inclusive differential cross section, from collisions of a hadronic projectile hp on 
a nucleus A, for the production of a hadron h that carries a longitudinal momentum fraction 
z, integrated over its transverse momentum. The differential leading particle multiplicity is 
then defined as 

rfAry'^^(^) _ 1 da'^'^'iz) 



dz ~(Tt;\A) dz ' 

where (T^"^'(A) is the inelastic cross section for hp — A scattering. In the limit of interest, 
zdN^^^ [z) / dz takes on the asymptotic form: 

/ ^ , ^ ^ = E / dx-D^/\-,Q')U„^{x,Q\ (6) 



involving the convolution of fa/hp{x, Q"^), the distribution of parton a in the hadron projectile 
hp, and D^^'^{z, Q^), the fragmentation function for parton a into hadron h. Eq. (|^) leads to 
a steep decrease of the hadron spectrum at large z because at ~ 1 GeV^ the structure 
functions drop as xfq/hp{x, Q"^) ~ (1 — a;)" with n ~ 2 — 3 and the dominant fragmentation 
functions zD'^^'^(z, Q"^) ~ (1 — 2) based on quark counting rules or at most are constant in low 
virtuality approximations ||T8|Jl9t| . The hard interactions that transfer transverse momenta 
to the fast partons are at low virtuality. We set the virtuality at = 1 GeV^. 

For all our calculations, we used the parton distributions of GRV |^0| since they are 
determined down to the low- virtuality range ~ 1 GeV^. Two types of fragmentation 
functions were used. One was determined by the EMC collaboration at < >= 
20 GeV^ and the other was calculated [|18| for low virtuality from the quark-gluon string 
model (QGSM) [|18[. The EMC fragmentation functions set a lower bound on our multiplicity 
estimates, since scaling violations will enhance fragmentation at low virtuality and large z. 
As a cross check, we fit the EMC data for to the higher twist form calculated in 

T^. We then evolved this down to = 1 GeV^. In this region in particular the Berger 



form also becomes nonzero at z = 1 which is consistent with the QGSM form. Numerically 
we found that the fitted form of was a factor of one to two times bigger than the QGSM 
form. We consider this discrepancy acceptable. 

All three fragmentation functions suppress light quark fragmentation into protons. The 
EMC data indicate a factor of 4 to 8 depletion of protons versus 7r+ in the largest multiplicity 



region of leading particles z ~ 0.2 — 0.3. The higher twist mechanism in |]T9| is inapplicable 
for proton fragmentation. The QGSM results [0 give zero fragmentation into protons. 

The limiting curve for proton (solid) and pion (dashed) projectiles are shown in Fig. 1 
using the QGSM fragmentation functions. One typically finds the leading behavior as z —>■ 1 
to be zdN^^^^{z)/dz ~ (1 — z)°''^''^p . Here ahjhp will be referred to as the leading exponent 
for produced species h from projectile hp with ax/hp denoting the exponent for the total 
leading particle multiplicity. For the proton projectile with QGSM fragmentation functions, 
we find a7r+/p ~ 4.4, 0;^-/^ ~ 5.3 and ax/p ~ 4.5. For the EMC fragmentation functions 
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axi'p ~ 6 for all species X except for the proton, where Op/p ~ 7, and a^/p ~ 6. For the pion 
projectile with QGSM fragmentation functions, we find Q;7r+/7r+ ~ 1-6, a^,- j^.^ ~ 2.8 and 
Q^T/7r+ ~ 1-8; for the EMC fragmentation functions Q;T/7r+ ~ 3.0. In all cases, multiplicity 
distributions for p and vr" projectiles are obtained by applying charge conjugation to those 
for the p and 7r+ projectiles, respectively. 

Let us now turn to what we might expect for the leading particle multiplicities. In a large 
nucleus, say A = 200, the relative transverse momentum that two leading partons acquire 



can be estimated following p2| through the px - broadening of the Drell-Yan spectrum. 
Using data fl^ for A(p|.)pA~2oo ~ 0.114 GeV^ and the average number of struck nucleons 
in the Drell-Yan events n(A=200) = 5.3 p3| : 



PTn A ~ 2^A{p^t)p AfiA/ifiA - l)U~2oo ~ 0.75 GcV/c. (7) 

This number should be considered as a lower limit since in QCD p'^ - broadening is larger 
for the interaction of partons with smaller virtualities |^ and in our case the virtuality 
is much smaller than in the Drell-Yan processes studied in |T^ for (M^+^_) = 30 GeV^. 
Within the z > 0.2 region, we estimate the coalescence probability into a pion, P^{p'^), for 
a qq pair with relative momentum p^^. It is easy to show in the constituent quark model 
that 

P:iPTj = \FMj\' (8) 

where Ft^{p'^) is the pion form factor. Using the Vector-Dominance-Model fit for F-j^lp"^), one 
finds for Eq. (|^) P^{0.5) ~ 0.25. Further suppression comes from the presence of gluons in 
the pion wave function which on average carry half of the pion momentum. 

For large A, the above mechanism for leading particle production goes as A"^/^. One 
finds a similar A-dependence for the fragmentation of a diquark to a nucleon. Another 



mechanism for the production of leading particles which is similar to |2^ is the filtering 
of color singlet small transverse size clusters of the projectile. This is analogous to the 
propagation of ultrarelativistic positronium through say a piece of lead. Using equations 



derived in |26| , one finds for large A tliat tliis meclianism leads to a multiplicity of leading 
pions (protons) oc A^'^^^ (A~^^^) . 

Another "background" to our mechanism is the coherent diffractive production of leading 
particles. Its contribution to zdN^^"^ /dz oc A^^^^ for A ^ oo and somewhat slower for 
A = 12 - 200 [0. 

Although we see no expectation for the asymptotic form in Eq. (|]) even for A = 200, we 
feel large nuclei may still inhibit coalescence. To test this claim, since the leading correction 
to the asymptotic behavior is due to diffractive events, it is preferable to perform data 
analyses which would explicitly exclude diffractive events by selecting the events where the 
nucleus breaks up. For this these asymptotic expectations can be compared with data from 
central hadron-nucleus collisions from Fermilab experiment E597 p8|Pj29||. 

The data were obtained using the Fermilab 30-inch hybrid bubble chamber spectrometer 
with associated downstream particle identifiers (DPI). The bubble chamber, in a 2T mag- 
netic field, provided a visual target and vertex detector as well as a spectrometer for the 
slower produced particles. The faster particles were momentum analyzed using the fringe 
field of the bubble chamber magnet and seven planes of proportional wire chambers and 
three drift chambers. In addition to being filled with liquid hydrogen, the bubble chamber 
also contained thin nuclear foils of Mg, Ag and Au. The beams consisted of 100 GeV/c 
p,p,7i~^ and Tr~ . Mass identification for these particles was provided by Cerenkov counters 
in the beam-line. Additional experimental details are given elsewhere . 



To make the relevant comparisons, leading particle spectra have been obtained by com- 
bining p/p and combining vr^ projectile data. To study the nuclear thickness dependence 
of the data, use is made of the known correlation between the number of nuclear collisions 
and the number of observed slow protons I^Oj. The number (up) of slow protons (with 
momentum less than 1.3 GeV/c) is determined by performing a visual ionization scan of 
each interaction occurring in the nuclear foil targets. The selection of > 1 automatically 
removes all diffractive events. 

In Fig. 2a the differential multiplicity zdN^"^ {z) / dz is shown for the data from tt~ A 
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/i+ (combined with n'^A h~) with = 1 or 2. In these figures, h~^{h~) indicates a 
positively (negatively) charged hadron, mostly pions. The dashed curve is a fit to the form 
(1 — 2;)" and demonstrates that this form is an acceptable representation of the data with 
a = 3.27 ± 0.48. Figs. 2b and 2c show the resulting values of a as a function of Up for the vr 
projectile to tt^/tt^ and p/p, respectively, with z > 0.2 and for produced hadrons having the 
same or opposite charge to that of the beam projectile. The solid (dashed) horizontal lines 
at the right side of the figure are the theoretical predictions using the QGSM fragmentation 
functions in Fig. 2b and the EMC fragmentation functions in Fig. 2c for the opposite (like) 
charge leading particle exponent. 

We next compare various integrated multiplicities from the experiment to the predictions. 
Let Ih/hp{zm) denote the integrated multiplicity of produced hadrons of type h with z > Zm 
from a projectile hp with h replaced by T when all the produced hadrons are summed. For 
the predictions, these are understood to be for A — 00. In table 1 the limiting predictions 
for the integrated multiplicities are compared with the experimental results for Up > 15. 
Fig. 2d shows the E597 data for the integrated multiplicity with Zm = 0.2 as a function 
of Up. The solid (dashed) horizontal lines at the right are the predictions from the QGSM 
fragmentation functions for p (it) projectiles. |^ The data seem to be consistent with the 
idea of limiting multiplicity with the magnitude as determined by the model. 

It is also interesting to examine multiplicity ratios. Denote the ratio as. 

For the ratio between protons and pions, only the EMC fragmentation functions are avail- 
able. The calculated values are (0.2) = 0.18 and (0.3) = 0.17. To 
examine the depletion of leading protons from a proton projectile, the lower bound estimates 



^ Using QGSM versus GRV proton distribution functions [31|, the predicted multiplicity is about 
the same. All the leading exponents are about two smaller, which is consistent with the slower 
decrease of the former as x — > 1. 
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from the EMC fragmentation functions are /p/p(0.2) = 0.009 and /p/p(0.3) = 0.002. The 
statistics on experimental data for p ^ p are insufficient to quote. It is worth noting that 
the forward multiphcity for large Hp is substantially larger than in eikonal models where 



usually it is assumed that the energy is equally divided between all exchanges, see e.g. [10 



In these models one should expect that almost no particles are produced for z > 0.2 for the 
case in which more than 5 nucleons have been struck. 

In our estimate of px - broadening we have used experimental data on Drell -Yan pro- 
duction at 400 GeV/c. No estimates of the energy dependence of the pt - broadening are 
available at the moment. The model used in the analysis 0] leads to energy independent 
energy losses. However this model assumes that the soft cut-off does not depend on energy. 
It seems natural to expect a certain increase of hardness of the soft interaction with energy 
leading to an increase of pt - broadening with an increase of energy. So fewer leading par- 
ticles will be produced in the central collisions as the energy of hA{AA) collisions increases 
from the SPS energy range to the RHIC or LHC energy range both because of further 
suppression of the coalescence effect and the scaling violations for the quark fragmentation 
functions. 

An implication of the current analysis is that in high energy central AA collisions very 
few baryons should be left in the fragmentation region. In fact our analysis indicates that the 
probability for a baryon to carry more than z = 0.2 of the projectile's momentum fraction is 
less than 1%. To conserve baryon number, a guess is that nucleons should move on average 
at least 4 units of rapidity to the central rapidity region for RHIC and beyond. 
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Foundation. 
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FIGURE CAPTIONS 

Figure 1: Limiting curves for leading n"^ production from p - A (solid) and tt"^ -A (dashed) 
collisions. 

Figure 2: Experimental results ior h — A collisions at lOOGeV/c: (a) Differential multi- 
plicity for Tr~A TT"^ combined with tt'^A tt" for events with Up — 1,2 {rip is the number 
of slow protons). The dashed curve is a fit to the form (1 — ^)"; (b) Leading exponent a for 
-kA — > TT, solid circles (open boxes) when the produced tt has the same (opposite) charge to 
that of the beam projectile; the horizontal lines on the right are the theoretical limits from 
QGSM fragmentation functions to the same (dashed) and opposite (solid) charge leading 
particle; (c) Leading exponent a for nA — > p(p), the horizontal lines on the right are the 
theoretical results from EMC fragmentation functions to hke (dashed) and opposite (solid) 
charge leading particle; (d) The integrated multiplicity of hadrons with > 0.2 for tt^ 
(solid circles) and proton (open boxes) beam projectiles as a function of Up. The horizontal 
lines on the right of the figure are the theoretical asymptotic limits for tt (dashed) and proton 
(solid) projectiles. 
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Table 1: The integrated multiplicities, Ifi/fi^[zm)^ of produced hadrons of type h with z > 
from a projectile hp. The first two columns give the predictions for the limiting case A — )• oo 
for the EMC and QGSM fragmentation functions, respectively, and the third column gives the 
experimental value for Up > 15, each for z^ = 0.2. The last three columns correspond to the 
case with z^ = 0.3. 



